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Abstract 

The numerical simulation of an air spark gap has been carried within two theoretical 
models. The kinetic one allowed us to calculate time dependencies for the residual 
resistance (0.2 - 0.4 Ohm for our selection of a circuit parameters), the spark gap channel 
width, the electron number density, the mobility, the conductivity, the ionization degree, 
the magnetic field in the discharge channel, the channel inductance and the electron drift 
velocity. Simulating a real circuit and taking into account a spark gap residual resistance 
demonstrates good agreement of both models with the experimental data, while that 
without taking into account this resistance overestimates the maximal current in the circuit 
by approximately 5%. 


1 Introduction 

The high-voltage spark gap switches are essential elements of the circuits containing magneto- 
cumulative generators, capacitive and inductive energy storages, electro-explosive circuit break¬ 
ers etc. Very frequently spark gaps affect sufficiently the operation of the circuits and one cannot 
satisfactorily describe main elements functioning in the absence of a quantitative model of the 
processes occurring in the spark gaps. This article is devoted to the comparison of two such 
models considering a real circuit, where the influence of a high-voltage spark gap results in a 
decrease of the maximal current for about 5%. 

We review in detail the model based on the time dependent method for the electron 
number density simulation in the air spark gap channel using an anisotropic solution of the 
kinetic equation [^. The solution was obtained taking into account vibrational, dissociation, 
excitation, and ionization losses and is applicable at held strengths E/po lying between 40 
and 450 V/cm-mm Hg. This model was successfully implemented for the experimental data 
interpretation in studying nanosecond high-voltage spark gaps at the upper region of the 
the specihed range of E/pq. 

We, on the contrary, apply the model Q in the lower region of the allowable held strengths 
E/po appearing in the spark gap at a voltage of about 15 -20 kV (that corresponds about 50 
V/cm-mm Hg). Apart from our main objective to describe experimentally observed inhuence 
of the spark gap residual resistance on the functioning of the whole circuit, we simulate vari¬ 
ous physical parameters of the spark gap channel and compare simulation results with those 
obtained using known model j^. 
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2 Physical models 

In simulation of all physical processes in a spark gap (dissociation, molecules ionization and 
excitation, electron recombination etc.) is reduced to the calculation of the electron number 
density depending on the voltage Us{t) and the field strength E{t) = Us{t)/d, where d is the 
gap length. We assume the field to be uniform along the discharge channel following j^, where 
it has been demonstrated that during a short enough time period (comparing to the time of 
a spark gap breakdown) the field distribution becomes nearly uniform. Also in [^, based on 
the numerical simulations of the discharge dynamics in the streamer model within 2D and ID 
approximations, it is stated that one can describe the processes in a spark gap channel with a 
reasonable accuracy within OD approximation. 

The electron number density is defined by the kinetic model on the assumption of an 
anisotropic distribution function and can be written as 


ne{t) = ne(0) exp^^ 

where no{t) = lcm~^ is the electron initial number density in the air and Q, 
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is the mean growth rate for the free electron number density. Following further z = 

1/2 / \5/6 

Va-As = 1.19 ■ ( & ) J/s, p is the air pressure in a spark gap, po is the 
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air pressure at sea level, is the excitation frequency of the air particles, Ae is the excitation 
energy, Si = 2.35 ■ 10“^® J is the ionization energy. ,^o = 2^1 which electrons are 

gaining energy in the applied electric field E, 


Vc = l-lln^/^cTc (—^ = 6.93 ■ 10^® f s ^ (3) 

\mej \PoJ 

is the elastic collision frequency of a free electron with the air molecules, e is the electron charge, 
rrie is the electron mass, is the number density of gas particles, cTc is the cross section for the 
inelastic collision of an electron with an air particle and it is gained from the experiment and 
depends on the air pressure [^. 

Knowing the time dependent electron number density, one can easily find the spark gap 
resistance 


m 


or the current flowing through it 


d 

esz(t)ne(t)pe(t)' 


(4) 


I{t) = a{t)E{t)s^{t) = esz{t)ne{t)fie{t)^^^, (5) 

where Pe{t) = is the electron mobility, Sz{t) = vrr^(t) is the current conductive area, 
a{t) = e/j,eit)ne{t) is the conductivity and r^(t) is the discharge channel radius. In most of 
papers on theoretical calculations one chooses a discharge channel width to be a constant 
value of about 10 microns [^. This agrees well with our own simulations. Let us give some 
additional reasons for our choice of a discharge channel model and its width; 1) if the 
width of a discharge channel is much greater than its radius, one can consider all the processes 
in the interelectrode gap to be one-dimensional. 2) The real value of the discharge channel 
radius defining the channel resistance, generally does not match the value of the luminous 
column radius which is registered experimentally. 3) The discharge channel radius should be 
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defined by the shortest distance from the channel axis to the space point where the plasma 
specihc conductivity is e = 2.7 times less than the maximal one. 4) One can consider the 
discharge channel width to be approximately constant [^, this statement was done on the basis 
of the numerical simulations in the streamer model within 2D approximation. We give below 
comparative simulation results for some physical quantities in the discharge channel calculated 
both at its constant and variable width and show that in the model Q the channel width 
depends on time, however for to estimate the residual resistance it is quite acceptable to choose 
it to be constant. 

To hnd a suppositional form of the time dependent width of the current conductive region 
we use an approximate expression 


r,(i) = \/2D(t)t, 


( 6 ) 


where D{t) is the electron diffusion coefficient and one can calculate it by two different formulas 
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and 1[^ 
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D{t) = (8) 

where Ut = ^ is the electron average thermal velocity and Umoi = 2.51 ■ 10^® cm“^ is the 
number density of the air molecules. We set an initial voltage on the spark gap to be 0.0001 V 
(if one sets it equal to zero, than there appears a singularity in the system of equations at an 
initial time moment due to a division by zero), the cathode-anode gap d = 1 cm, the pressure 
p = 120 kPa and the cross section cTc = 0.2 ■ 10“^^ m^. There is the only empirical parameter 
in the model described - the elastic collision frequency of a free electron for which we use the 
experimental data printed in [^. 

In the second model we use, the air gap is replaced by a current generator and the spark 
gap is described by a nonlinear second order differential equation of 


y{t)y{t) + y{t){l - y{t)) + 6(1 - y{t) - y{t))ip{ay{t)) = 0, 

yit) = Vi.ay{t)) = a{E/pf, (9) 

which is solved together with the circuit equation. Uq is the breakdown voltage of the spark 
gap, a is the impact ionization coefficient, we use in our calculations b = 24.472, a = 15.4 ■ 10® 
the other coefficients can be found in [^. Unlike Q it is impossible to get time dependencies 
for the electron number density ne{t), the discharge channel width rz{t) etc. in model |^. 

3 Simulation results 

To compare the predictive capabilities of the models described, we simulated a spark gap as a 
part of a simple electrical circuit consisting of 19.2 mcF capacitor, 300 mcHn inductance, 0.15 
Ohm resistance and initial 15 kV voltage on the capacitor. We compare simulation results with 
the experimental data on the current in a real circuit with the same parameters. 

There are four current graphs on Fig.l ; the experimental, the simulated by model 0.by 
model and under the assumption of zero spark gap resistance. The last graph demonstrates 
the importance of taking into account a spark gap resistance because there is an overestimation 
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Figure 1: Currents: experimental (solid, blue), simulated by model Q (thick,red), by model 
(dashed,green), under the assumption of zero spark gap resistance (dotted,yellow). 


of the maximal current in the circuit by approximately 5% without it. The currents simulated 
by models are almost identical with the experimental results in the hrst oscillation half¬ 
period in the system which exceeds 200 microseconds. 
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Figure 2: Resistances: simulated by models (thick) and |[^| (dashed) on the time intervals 
(left) from 0 to 10 /iS and (right) from 18 to 200 ps. 


The only drawback is the 8 microsecond current shift in the model with respect to the 
experimental data. This circumstance is explained by the fact that the model |[^| originally is not 
kinetic, but as is known, kinetic processes dominate in the spark gap in the hrst microsecond. 
One can see from Fig.2 that the resistance subsides slower in than in Q and tends to decrease 
to a small nonzero value in the model 1|^|. In the model Q the resistance falls to 0.4-0.2 Ohm 
(depending on the circuit parameters) and then begins to increase, which may indicate spark 
gap locking. Since we have the ability to compare with a real experiment, we conclude from 
the above that the model Q is closer to the real situation than |^, although they both predict 
very close values of the resistance and the current. 

The kinetic nature of the discharge process at the hrst microsecond represents itself an 
essential advantage of the model Q with respect to [^. The model also allows one to calculate 
numerically the evolution of various physical characteristics during the discharge process. Fig.3 
presents the simulated graphs for the spark gap channel width r^(f), the electron number 
density ne{t), the mobility the conductivity cr(f), the collision frequency Uc, the ionization 
degree ne{t)/nmoi, the magnetic held H{t) in the discharge channel, the channel inductance L{t) 
and the electron drift velocity u}e{t) = fie{t)E{t). Solid lines correspond to simulations with 
f'ziit) = where D{t) is determined by Eq. (7), dotted ones to that with = const = 10 
microns (we simulated also larger values of r^, however they almost did not ahect the time 
dependencies of the other physical quantities) and dashed with rz 2 {t) = Tz{t), where D{t) is 
determined by Eq. (8). 

As seen from Fig.3 the electron mobility, the elastic collision frequency of a free electron 
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Figure 3; From left to right and from top to bottom: time dependencies for the spark gap 
channel width rz{t), the electron number density ne{t), the mobility the conductivity 

cr(f), the collision frequency the ionization degree nf.{t)/nmoh the magnetic field H[t), the 
channel inductance L{t) and the electron drift velocity simulated by the model [^. 


with the air molecules, the channel inductance, the electron drift velocity as well as the ohmic 
resistance are almost independent of r^(f), although the electron density and conductivity are 
very sensitive to it. Namely the graphs ne{t) and a{t) show that at = const the density 
and the conductivity reach their maxima at the current maximum, while for rz{t) even slightly 
varying with respect to = const = 10 microns their maxima shift in time to the beginning 
of the process, that is more consistent with the actual physical picture of the process. 

4 Conclusion 

We have simulated the residual ohmic resistance in the channel of an air spark gap using 
two different physical models and estimated its value as 0.2 - 0.4 Ohm in the case of the 
simple electrical circuit. The simulated values agree well between themselves and with the 
experimental data. We also show that the time dependent width of the discharge channel 
allows one to describe the physical processes in a spark gap much more correctly than the 
constant one. 
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